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Abstract 

Supposing only Lorentz and the gauge invariances of the Lagrangian, 
we derive energy and angular distributions for n^ii^ii^ lepton 

flavor violating decay process. Using these results, we discuss methods 
to determine the parameters associated with the lepton flavor violating 
interactions. 



* akihiro® eken . phy s . nagoy a-u . ac . j p 
^ sanda@ kanagawa- u . ac . j p 



1 



1 Introduction 



Lepton flavor violating (LFV) decays, for example, /j, 67, r ^7, fi-e 
conversion and t ^ Sfj, are being studied extensively by experimentists [l]- 
[1], and by theorists [5]- [8]. Especially, a large number of r decay events 
are collected in B Factories and we hope that LFV decay mode of r may be 
found in the near future. In that case, the Super B Factory 9j might give 
us the large number of r — > 3;U events. We might also hope to observe the 
energy and polarization distributions. 

In r — > 3/i decay, we can measure all the energies and directions of the 
final state muons. For definiteness, let's say that we want to investigate the 
r"*" decay. B Factories generate back to back t~^t~ pairs. As mentioned in 
Section El the polarization of can be observed statistically by correlation 
of momenta of decay products of both r"*" and r~ [lOj. For fi^ e'^e~^e~ 
process, the differential branching ratio which is a function of energies of 
2 positrons in the final state and the polarization of in initial state is 
already derived in Ref. [5]. The structure of weak interaction which causes 
the /i — > evu decay was investigated by Michel [11]. He introduced Michel 
parameters which proved to be very useful. We follow the similar strategies 
and formulate the method with which we can probe the structure of LFV 
interactions starting from a general Lagrangian. 

The general Lagrangian for decaying to fj,~^ fj,~ supposing Lorentz 
invar iance and gauge invar iance is written as: 



where and rrir are the masses of the fj. and t^, respectively, Gp is 
the Fermi constant. {tl^P'L-, I^r\ and {fR,flR,fj,L} are the Dirac spinors 
with the helicity operators, (1 it 75)/2, respectively. For example, fpt = 
f(l-75)/2 and fiR = (l + 75)/i/2. a"^ = |(7"7^-7^7"), D° = 9° + ie^°, 
p^afB _ Qaj^i3 _ Qli j^a ^ j^a ^-j^e photon field and e = — |e| is the electron 
charge. C-y represents r — > /i7 interaction as well as the terms present in 



C = Cpp -\- C. 
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Figure 1: These are the diagrams produced from L. 




Figure 2: These are the diagrams produced from Lpp. 

ordinary QED. These terms generate diagrams in Fig. [TJ So, Ai^ and A^^ 
are the coefhcients of interactions in which the intermediate photon has the 
left polarization, and the right polarization, respectively. The coefficients 
Q\ ~ 56 in C.FF are the coefficients of various 4- Fermi type interactions. 
These terms generate diagrams in Fig. [2l Using the Fierz transformation, 
it is shown in the Appendix A that Eq. ([T]) is the general Lagrangian. 

For convenience, we summarize our major results, here assuming that 
enough events are collected. 

1. We have shown that the absolute values of the coefficients |(7ip/16 + 
|53p + b2|Vl6 + b4p, l^sP + beP, + Re.\g^Al\^Re\g^A*j^ 
and Re\g^A*j^ + Re\g^A*^ can be obtained by measuring the energy 
distributions of decay products of r^. 

2. We have shown that the absolute values of the coefficients |(7ip/16 + 
bap, I52IVI6+I54P, bsl, l^el, \Ar\, I^lI, -Re[c/3^I], Re\g^A*^, -Re[55^I], 
i?e[fif6^y , lm{gj,A*^^lm\g^A*j^ and lm\g^A*^^Im\gf^A*j^ can be ob- 
tained by measuring the momentum distribution of decay products of 
T~ in addition to the energy distributions of decay products of r+ . 

3. In some suitable cases, we can determine l^iil, \g2\-, l^al and \gi\^ inde- 
pendently, and phases wgXg^A*^^, &Tg\g'>,A*j}, w:g\gi^A*j^ and wg\g^A*{\ 
as shown in section [9l 

4. It is interesting to note that there is a possibility that some information 
on r — > /i7 could be obtained before r — > /i7 decay is measured. For 
example, i?e[54A|j] and i?e[53^2] arise from the interference terms 
between the four Fermi interaction amplitudes and the r — > |U7 decay 
amplitudes. The decay rate for r ^ //7 is quadratic in \Ar\ or \Ai\ 
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while the observables given in Re[g4^A*fi\ and Re[g3A*j^] are hnear in 
\Afi\ and So the interference effect may be seen before the decay 

rate for r ^ /U7 is seen. 

This paper is organized as follows. Section [2] gives the differential branch- 
ing ratio of r"*" — > decay. Also, a general formula which yields in- 
formation on the current structure for t+ decay is derived. In 
section [3l we pick energetic one of the two n the final state and analyze 
its energy dependence. In section HI we discuss physics implications of our 
results using only the results of section [3l In section O we give the energy 
dependence of fi~ . In section [6l we discuss physics implications of our re- 
sults using only the results of section \5[ In section \7\ we discuss physics 
implications of our results using the results of sections U] and [6l In section 
IHl we give polarization dependence of the branching ratio in addition to 
the energetic /u+ energy dependence and analyze the results. In section [H 
we discuss physics implications of our results using the results derived until 
previous section. Section [TOl contains the concluding remarks. 



2 General Formula 

In this section, we derive the differential branching ratio for t+ jj,^ 
decay including r+ polarization, and the general formula for the observables 
which are relevant for the actual experimental situation. 

The final state contains two mesons. The one with higher energy is 
denoted as /ii. The other is denoted as /i2- M~ is denoted as 

In the real experiment, the processes we want to detect are 

e^e~ T~^{s^) T~{s~) —^i'T- + a + anything 

^ Ml M2 M3 

where a is a particle which has the charge —1. So we must calculate the 
differential cross section for these processes. 

Here, we use the center of mass frame of e~^e~ initial state, which we call 
frame 1. In the rest frame of t~^, which we name frame 3, the momenta of ^ui, 
IJ,2, M3 and t+ are denoted as pi = {Ei,pi), p2 = P2), P3 = {E3, pa) and 
Pr = {mr, 0), and finally the polarization of r"*" is denoted as denotes 
the polarization of r~ in the rest frame of r^. Here, s'^ -pr = s~ -kr = and 
(s^)^ = — 1 where k-j- is the momentum of t~ in the t~ rest frame which we 
call frame 2. 

For the definiteness, we set the relations between frame 1, 2, 3, as follows 
and they are depicted in Fig. [3l 

Frame 1 is the center of mass frame of e~^e~ . Defining the t+ momentum 
as p!^ and the initial state positron momentum as p'^+, we set the z 
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Figure 3: The relation between frame 1, frame 2 and frame 3. rj is the angle 
between and p'^. 

direction in this frame as the direction of Pg+, y direction is the same 
as X p'^ and the x direction as {p'^+ x p!^) x p'^+. 

Frame 2 is the t~ rest frame. We set the z direction in this frame as the 
direction of p'^. The y direction is the same as that of the frame 1. 
The X direction id defined by (p^^ x p'^) x p'^. 

Frame 3 is the r"*" rest frame, We set the directions in this frame are the 
same as that of the frame 2. 

There are four ingredients that are used to compute the differential cross 
section for the process shown in Eq. (jlj). 

1. The narrow width approximation where we approximate 

e'^e^ T^(s^) r^(s^) ^ + o- + anything 

^ Ml M2 M3 

is described in appendix IB. 3i 

2. The differential cross section for e^e~ j* ^ t^(s^)t^{s^) is pre- 
sented in appendix lB.il 
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3. The differential branching ratio for r (s ) — > i^,- + o + anything is 
presented in appendix IB .21 

4. The differential branching ratio for t^(s^) ^1/^2/^3 is given below. 
The differential cross section for this process da is written as [12] 



da 



dQ dxidx2dil.T-dip d^ka 

da(^e^e^ r+(s+)T^(s")) 



S E 



dn 



±s+,±s- 

dBr{T^ {s~) Vt + a + anything) 



(5) 



X 



d^ka 

dBr(T+{s-^) ^i/U2/i3) 
^ i_ 

dxidx2dQrdip 

S implies sum over polarizations, ka is the momentum of the particle a in 
T~ rest frame. The definitions of ^1^ and '0 are written in next paragraph. 
We have computed the branching ratio for r+(s+) A*i//2A*3 including 
dependences. However, for now, we confine our discussion where we can 
approximate m^/rrij- = 0. The result is 

(i.Br(r+(j;+) ^^21^-3) 
dxidx2dil.T-dip 

3 _ r (6) 

= -^Brir ^ Uvu) ^0(2:1, X2) + 2_, ' Pj^Kxi, X2) , 

i 

where the definitions of Go(xi,X2) and G|(xi,X2) are written in Appendix 
ICl Br(r — > ^uu) is the branching ratio of r"^ — > ^'^v^Dt decay. 

2Ei 2E2 2Es , . 

xi = , X2 = , X3 = . (7) 

nir TTlr rrij- 

Note that xi, X2 and 2:3 take values between and 1, and xi + X2 + x^ = 2. 
Vl-T is the solid angle defined in Fig. [Hand ip is the angle between the P2-P3 
plane and z-pa plane as defined in Fig. El Pj = {pi, p2, Pi x P2} where 
Pi = Pi/|Pi| and P2=P2/|P2|- 
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frame 3 
z P3 




Figure 4: Definition of Q and (/> in frame 3 




Figure 5: V is the angle between P2-P3 plane and P3-2 plane. ^ is the 
angle between p3 and pi. x is the angle between p3 and p2. < ^ < vr. 
< X ^ 27r. 
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Substituting the concrete representations, Eq. ([5]) becomes 
da 



dVL dxidx2dO,rdilj d^ka 
iq 



2^ A«2 



2 sm" r/ 
( 1 + cos 7? H TT- 

T 



2 ■ 1 

/-, 2 sm^??, 
+ (l + cos T] —)s:,s. 
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2 )"z 



— sin^ T] SySy + (1 H — 2 ) ^ ■s^ -s^ 



sin 2?7 

7 



■^2 '^a; "I" '5 2 



— Sr(r ^ //z^i>) Go(a;i,X2) + • PiG|(xi,X2^ 



Br(T -^Ur + a + anything) — \— \G1{ya) - s • k^Gf (y, 



(8) 



Carrying out the summation for ±5+ and its using s^s^ = 5ij, Eq. 
becomes 



da 



d^l dxidx2d^lrdip d^ka 



4i3r(r — > ^i'v)Br[T i/^- + a + anything) — ^r-^ — ^— 



2 

Go{xi,X2)Gl{ya){l + cos^ r, + ^^) 

. 2 

- ^ X2)G5(2/a){ (1 + cos^ r/ - ^^) fc^^P^^ - (5'^ siv? r, kayPi 



sin 2r] 



+ (1 + 2) ^ kaxPix ikaxPiz + kazPi 

(9) 

This expression is the general formula of the process (jH) as long as we stay 
sufficiently away from singularity at = 0. In following sections, we'll 
start analyzing from this expression. 



3 ^0(3:1): Energy Dependence of fii 

In this section, we derive the formulae convenient for investigating the struc- 
ture of the LEV four Eermi interactions, using the ^1 energy dependence of 
the differential branching ratio for r — > 3/i decay mode. 

Our first priority is to discuss the observable which is easier to detect 
and analyze. So here, we integrate the polarization dependence as follows. 
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First, we integrate over dQdi^rdijj, 
da 

dxidx2 d^ka 

= Br{T iivv)BriT~ ^ Vr + a + anything)^ — ^GQ{xi,X2)G\{ya){2 + 

(10) 



Next, we integrate over d'^ka, 

da 
dxidx2 

32TTa^(3 . 1 , „ , -^r. / - , N (11) 

= K — (2 H — T^)Br(T ^ fjii/u)Br[T —> i/,- + a + anything) ^ ' 

X Go{xi,X2). 

Eq. ([TT]) allow us to obtain Go{xi, X2)- It has the arguments xi and X2- 
Even expression, however, is pretty complicated for a discussion here. We 
thus discuss only the xi dependence integrating over X2- For the physical 
region, we found that the effect of neglecting muon masses in the differential 
branching ratio for the r — > 3/i decay introduces an error of 0{2mfj_/m-r)- 
Now we define 

{\9l\ , I |2\ , /\92\ , I |2\ 

«+ = [-^ + \9^\ ) + + ) 



K = \95? + \9g? 



c+ = |eAfi|2 + le^Lp 

d+ = -{Re[g^eAl\+ Re[gieA\]) 

e+ = -{Re[gQeA\] + Re[g^eAl]). 

Then, Gq integrated over X2 is 

Gq{xi) = j dX2Goixi,X2) 

1 — X'l 

= ^Ji(2xi-l){6(l-xi)(2xi-l) 

+ /5ia(6xi -5)(2xi -1) 
+ ^/Oi6(3xi-2)(xi-l)} 



(12) 



^ 4^ f 2(2x1 - l)ixt - a;i + 1 



+ 3(2x?-2xi + l)log[-^]|, 

1 — Xl J 



3^1 (1-xi) ' ' ' ' ^n-xr 

(13) 



where 
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Ji = -(10a+ + 76+ + 72d+ + 54e+) 
3 

Jipia = ^(4a+ + 6+ + 48d+ + 12e+) (14) 
9 

-^i/^ib = +8e+). 

Wc have defined pia and to take values between and 1. We now discuss 
how c+, Ji, pia and pi^ can be determined from the xi dependence of Gq. 
First, to determine c+, it is convenient to define the function 

Fi{xi) = l{l-Xi)Go{xi). (15) 

o 

By choosing the kinematics such that xi ^ 1 for the function Fi{xi), we 
can obtain c+. as 

Fi{xi)l^^^ = c+(2xi - l){xl - XI + 1)1^^^^ = c+. (16) 

Next, we subtract the term containing the coefficient c+ and define an- 
other function, 

Go{xi) — (c+,xi term) 



F2{xi) 



(2x1 - 1) 

= ^Ji{6(l-xi)(2a;i-l) ^^^^ 
+ pi„(6xi-5)(2xi-l) 
+ ^pi6(3xi-2)(xi-l)}, 

where 

(c+,xi term) = 4^^| 2(2xi -^l)(xf - xi + 1) ^ ^^^^^ _ + l)log[3^] 

(18)' 

We can then determine the parameters Ji, pia and pib from the shape of 
^2(3^1) using 

F2{xi)dxi = — Ji, (19) 

1/2 

i^2(xi)|^^^, = ^Jl/9i„ (20) 



and 



respectively. 



F2ixi)l^^i=PiPib, (21) 



" 2 
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4 Physics from Go{xi) Distribution of fii 



The conclusion of previous section is how to determine the parameters, c+, 
"^1) Pia, Pib- Using only these parameters, o+, 6+, d+ and e+ can't be 
determined. However, if we take some special cases which are explained 
bellow, we can restrict the allowed regions of a+, 6+, c/+ and e+. 



4.1 What happens if one of a_|_, c+, d^, e+ =0 

If wc impose that one of a+, b^, c+, d+ and e+ is zero from other experi- 
mental results or supposing some specific models, the results are as follows. 
When a+ = 0, also d+ = and 

Ji = ^(76+ + 54e+) 

Jipia = ^{b+ + 12e+) (22) 
9 

-^iPi6 = -(6+ + 8e+). 
When 6+ = 0, also e+ = and 

Jl = ^(10a+ + 72d+) 

JiPia = ^(4a+ + 48d+) (23) 
Jipib = 0. 

When c+ = 0, also d+ = e+ = and 

Jl = ^(10a+ + 76+) 

JiPia = \{^a+ + b+) (24) 
9 

JiPib = -:b+. 



When d+ = 0, 
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Ji = ^(10a+ + 76+ + 54e+) 

Jipia = ^(4a+ + 6+ + 12e+) (25) 
9 

-^iPi& = -(6+ + 8e+). 
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When e+ = 0, 

Ji = \{10a+ + 7b+ + 72d+) 

^ipia = ^(4a+ + 6+ + 48d+) 
9 

Jipib = -^K- 

If any one of a+, 6+, c+, (i+ and e+ is zero, we can determine a+, 6_| 
d+ and e+ from only the xi distribution. 

4.2 Results if either a+ 0, 6+ 7^ 0, (i+ 7^ or e+ 7^ 
Similarly, When only a+ 7^ 0, 

Ji = ya+ 



3 
5 

Pib = 0. 



When only 7^ 0, 



3 

27 

When only a+, c+, d+ 7^ but a+ is negligible, 

Ji = 24d+ 

Pife = 0. 

When only 6+, c+, e+ 7^ but 6+ is negligible, 

Ji = 18e+ 
1 

Pia-3 
P16 = 1- 

In these case, pia and pi;, are independent from a+, 6+, cZ+ and e+. 
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4.3 What happens ii t ^ jj^ dominates 

If r ^ fij type interaction which has coefficient c+ is much larger than other 
interactions, only c+, and e+ may be determined since (i+ and e+ are 
not quadratic but linier of 4-fermi type interactions. In this case, 

Ji = 6(4(i+ + 3e+) 
Jipia = 6(4d+ + e+) (31) 
Ji/9ib = 18e+. 

From Eqs. (l2M|) and (l236|) in appendix [Dl 

4 ^ 1 / Ji(i-pib) 



c+ c+ V 24 
^ ^ 1 fJmbV 



(32) 



So, we can determine a+ and 6+ lower limits, though we can't determine 
and 6+, directly. 



4.4 Lower hmits of c+ when r — is highly suppressed 

We give three types of c+ lower limits. Especially, these are very useful 
when r — > p'j is highly suppressed compared with r ^ 3^ and 0+ cannot 
determine directly. 

First, we give the relations, 







= 6+ + 8e+ 


3Ji - bJipia 


- 2JiPi6 


= -12(4(i+ + e_ 


3Ji — 3Jipia — 


22 


= 4(a+ - 2e+) 


3Ji + Jipia - 




= 12(a+ + 8d+) 


6Ji - lOJiPia - 


8 , 


= 3(6+ - 32d+) 


3Ji — 3Jipia 


- 2Jl/9lb 


= 4a+ + 6+ 



(33) 



from Eqs. (fTil) . 
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Using Eqs. ([SD, (l23i]) and (1236]) . 



1 (3Ji - 5Ji/)ia - 2 Jipib)2 (4d+ + e+) 



122 X 5 - 3Jipi„ - 2Jipib 5(4a+ + 6+) 

(4d+ + e+)2 



< 



5(4^ + ^ 



(34) 



(4|^ + 1 
5(4^ + 1 



<c+. 

This becomes equality when a+c^ — = 0, 6+0+ — e+ = and = e+. 

When a+ c+, if 6J1 — lOJipia — |Jipife = 3(6+ — 32ci+) < and no 
special cancelation between 6+ and 32(i+, then 6+ — 32(i+ = C'(— 32d+) and 
6+ < 0(32^+) <C 4a+. So, using Eq. (l2Mll . 



1 (6J1 - lOJlpla - I JlPlfc)^ 

256 X 32 3Ji - 3Jipia - 2Jipib 



32(i+)2 3224 



256(4a+ + 6+) " V256 x 4a 



(35) 



We note here that 4a+ + 6+ is similar to 4(a+ — 2e+) and 4(a+ + 8d+) in this 
situation. So, we have similar result if the denominator of left hand side of 
Eq. (f35]) . 3Ji — 3Ji/9ia — 2Jipib = 4a+ + 6+ is exchanged by 3Ji — SJipia — 
fJiPib = 4(a+ - 2e+) or (3Ji + Jipia - fJipib)/^ = 4(a + 8d+). 

Similarly, when 6+ ^ c+, if 3Ji — SJipia — ^JiPib = 4(a+ — 2e+) < 
and no special cancelation between a+ and 2e+, then a+ — 2e+ = C'(— 2e+) 
and a+ < 0(26+) <C 6+/4. So, using Eq. 



1 (3Ji - SJipia - "^Jipibf 

82 3Ji - 3Jipia - 2 Ji/9ib 



4(4a+ + 6+) - V 46 



(36) 



We note here that a+ + 6+/4 is similar to a+ — 2e+ and a+ + 8(i+ in this 
situation. So, we have similar result if the denominator of left hand side in 
Eq. ([36]) . 3Ji— 3Ji/9ia — 2Ji/9i5 = 4a++6+ is exchanged by ^Jipib = 6++8e+ 
or (6J1 - lOJipia - |JiPi6)/3 = 6+ - 32d+. 

If we suppose that 6+ = from some special models or other experimen- 
tal results, we have c+ lower limit from (f23]l and p34p . 

1 (3Ji - Uipia? _ 4 < (37) 



1728 Ji - Jipia a+ 
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Similarly, from (I22p and (12360 . if we suppose that a+ = from some 
special models or other experimental results, 

J_(3J,-14W = 4 < (38) 
1080 Ji - SJipia ^+ ~ 

4.5 A bound for d+ 

4.5.1 Without using information c_|_ 

We obtain a bound for d+ using information from Gq{xi) measurement 
without using information on C-|_. 
Prom the Eq. P^ . 

l_ 9 + 3pia-10pib ^ 8 , a+ , . 

36 3 - 3pia - 2/916 4a+ + + 4a+ + 6+ ' ^ ^ 

If d+ = 0, 

19 + 3pia - Wpib 1 . , 

< — — < - 40 

- 36 3 - 3pia - 2pib - 4 ^ ' 



since a+, 6+ > 0. So, if 



then < 0; if 



then (i+ > 0; and if 



then, from Eq. (|33 



< 0, (41) 

3 - Spia - 2pib 
9 + 3pia - lOpib 



3 - 3pia - 2pib 



> 9, (42) 



0< ^ + ^f""^;^" <9, (43) 
3 - 3/>ia - 2/>ib 



_ 3 - 3p,a - 2p,b < . < J 3 - 3p,a - 2p,b (4^^ 
32 - + - i 32 ^ ^ 

Furthermore, we can determine d+ allowed region as follows. Using the 

fact 

< < 7' (45) 

~ 4a+ + 6+ - 4' ^ ^ 

which become the equation when a+ = and 6+ = 0, respectively, 

8 ^ ^1 9 + 3pi, - lOpife ^ 8 1 

4a+ + 6+ + - 36 3 - 3/9ia - 2pife " 4a+ + 6+ + 4' ^ 

The solution of this inequality about d+ becomes 

-9Ji + ISJipia + AJipib ^ , ^ 9Ji + 3Jipia - lOJipib 

144 - + - 288 ■ K ) 

This becomes equality, when 6+ = and o+ = 0, respectively. 
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4.5.2 Using information c+ 

We obtain a bound for using information from Gq{xi) measurement with 
using information on c+. 
Prom Eq. ([HD, 

a+c+ + 8d+c+ = -^(9 + 3pia - Wpib)c+. (48) 
So 



This becomes 



"^^ (9 + 3pia - 10pib)c+ > 4 + 8d+c+ (49) 



36 

since a+c_|_ — > as explained in appendix [D1 Using c-f. > 
Ji ^ ^„ , fd 



(9 + 3/>ia - Wpib) > + 4^ - 16. (50) 



36c+ 

The solution of this inequality about d+ becomes 



-4c+ - c+ J ^(9 + 3pia - Wpib) + 16 



<d+< -4c+ + c+ J ^(9 + 3pia - Wpib) + 16. 



(51) 



So we have another limit of 

4.6 A bound for e+ 

4.6.1 Without using information c+ 

We obtain a bound for e+ using information from Go(xi) measurement 
without using information on 0+. 

Similar to the subsubsection S^TTl 

Pib ^ . I ^+ .t-o^ 



9 3 - 3pia - 2pib 4a+ + 6+ 4a+ + b. 
and if e+ = 0, 

0<1 ^1^ <l. (53) 

- 9 3 - 3pia - 2pu - ^ ' 

This becomes equation when 6+ = and a+ = respectively. So, if 

Pib 



3 - 3pia - 2pi6 



< 0, (54) 
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then e+ < 0; if 



3 - 3pia - 2pib 4 

then e+ > 0; and if 



> (55) 



0< , , <l (56) 



then, from Eq. (i33]l . 



3-3pia-2/)ife 3-3pia-2/)ib , . 

--'1 ^ < e+ < Ji . [01) 

Furthermore, we can determine e.^ allowed region as follows. Using the 
fact 

< , ^+ , < 1, (58) 
4o+ + 6+ ^ ^ 

which become the equation when = and a+ = 0, respectively, 

8 4 pifc ^ 8 

4a+ + 6+''+ - 93-3pia-2/)ife - 4a+ + 6+''+^ • ^ ^ 

The solution of this inequality about e+ becomes 

JlPlb _ 3Ji — 3Jipia — 2Jipih ^ ^ '/iPlfc /gg\ 

18 8 - ^ - 18 ■ 

This becomes equality, when 0+ = and 64. = 0, respectively. 

4.6.2 Using information c+ 

We obtain a bound for e+ using information from Go{xi) measurement with 
using information on c_|-. 

Similar to previous subsection, from Eq. (fUj) . 

4 

6+C+ + 8e+c+ = - JipibC+. (61) 

This becomes 

^JiPib+c+ >e\ + 8e+c+ (62) 
since 6+0+ — > as explained in appendix [Dl Using c+ > 0, 

^'^^^^^>f^ + 4V-16. (63) 



9 c+ \ c+ 
The solution of this inequality about e+ becomes 



-4c+ - 2c+ + 4 < e+ < -4c+ + 2c+ J ^ + 4. (64) 



So we have another limit of e+. 
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5 ^0(3:3): Energy Dependence of fi^ 



In the previous section, we have discussed how the values of parameters c+, 
•-^1) Pia and pib can be obtained from Go but we have not determined the 
values of a+, 64., d+ and e+. This will be the subject of this section. 



I3 = 2 - 211 




Figure 6: The physical region for t+ PiP2P^ is the area which is included 
in curbed line and besides the right side of the line X3 = 2 — 2xi. Near the 
line = 1, the differential branching ratio is singular since the intermediate 
state photon becomes real. 

In the previous section, we gave xi dependence for Gq. In this section, 
to determine the values of a+, (i+ and e+, we give x^ dependence of 
differential branching ratio. As in the previous section, we do not discuss 
the polarizations. For most of the analysis, muon mass dependence can 
be neglected. But, as we shall see, in some part of the phase space, 
dependence must be taken into account. 

To study the formula of the Gq as function of X3, we give the following 
prescription. We can get the X3 dependence of Go(xi, X2) by using xi-\-X2 + 
X3 = 2. From the condition xi > X2 and the relation xi + X2 + X3 = 2, 

^1 > 1 - y . (65) 

So, if all muon masses can be neglected, we have 1/2 < xi < 1 and < X3 < 
1 as shown in Fig. El The integration of xi in the X1-X3 plane is given by 

1 

j dxiGo{xi,X3). (66) 
^ 2 

When xi = 1 in massless limit of muons, the photon in Fig. [T] becomes 
on shell and the propagator becomes singular. This leads a divergence of 
xi = 1 in differential branching ratio for r ^ 3/i decay. 
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So we have to be more careful with the range of at xi including the muon 
mass dependence. The on shell constraint P2 = m'j^ leads to 



^ + l-xi-X3 + ^-^Pi-P3/^?^/^i^ = 0, (67) 



where 6 = 2m^/mr- The true integration range for xi is given by Eqs. ()65p 
and ()67p with — 1 < pi • Ps < 1 as shown in Fig. [6l We find it convenient 
to approximate the domain of integration as 1 — ^ < xi < 1 — (^)^. 
Here, the upper bound of xi is decided as follows. First, we calculate the 
total branching ratio of c+ sector integrated on the true physical region in 
0(5^). Next, we calculate the total branching ratio of c+ sector integrated 
on the approximated region. Finally, we set the approximated upper bound 
to match these two total branching ratios. This approximation allows us 
to avoid the divergent region. With this approximation, we give the X3 
dependence of Go as 

i-(f)2 

Go{xs) = j dxi Go(xi,X3) 

2 

if 2 1 

= 2'^3a;3S 3^3(1 - a^s) + P:iaXz{xz - -) + P'ibixz - 1)(2;3 - -) 

+ 4c+ 1 2x3 - 3x^ + 2(2x^ - 2x3 + l)log[ 



-3^/xi 



where 



J3 = ^(12a+ + 56+ + 144d+ + 36eH 



(68) 



(69) 



JzPia = 6+ + 12e+ 
J'iP'ib = 48(i+. 

First, to determine c+, it is convenient to define the function 

F3(X3) = (70) 



By choosing the kinematics such that X3 ^ for the function F^{x2), we 
can obtain c+ as 

F3(X3)L3^0 = C+. (71) 

Next, we define the function 

Co (3^3) - {c+,x-i term) 

If"' 2 1 1 (^2) 

2 "^3 { Xz{l - X3) +p3aX3(x3 - -) + /)36(x3 - l)(x3 " 3) 



^^4(2:3) 
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where 

(c+ , X3 term) = 4c+ 1 2x3 - + 2{2xl -2x3 + 1) log [^|^] | . (73) 

From the function F/^^x^), we can obtain the values of J3, p^a and ^35 by 
fitting the experimental data as follows: 



1 

dx-iFi{x'i) = r^h 







1 (74) 



F4(X3)U.3=1 = -J3P3a 
^4(2:3)1x3=0 = ^J-iPZb- 

In addition to Ji, pia and /Oi;, from previous section, the parameters a+, 
6+, dj^ and e+ can now be determined. For example, 

11 2 

«+ = ^«^3(l + gPSa - /?3fe) - - JlPlft 

4 

= - 2J3/33a 

3 (75) 

e+ = -^JsPsa - gJlPlb- 

Prom the differential cross section which are the functions of xi and X2, 
xi and X3 or X2 and X3, all we can determine are the quantities a+, 6+, c+, d+, e+. 
Now, we will discuss some specific cases in which we can go further. 



6 Physics from Go{x^) Distribution of fi^ 

Similar to the section HJ we have some bounds for a+, d^ and 64. from 
the analysis of previous section. 

6.1 What can we conclude if one of a+, 6+, c+, e+ =0 

If we impose that one of a^, 6+, c+, d^ and e+ is zero from other experi- 
mental results or supposing some specific models, the results are as follows. 
When a+ = 0, also d+ = and 

J3 = ^(56+ + 36e+) 

J3P3a = b+ + 12e+ (76) 
JsPsb = 0. 
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When 6+ = 0, also e+ = and 



J3 = 4(a+ + 12d+) 

J3P3a = (77) 



J3P3b = 48(i+. 
When c+ = 0, also d+ = e+ = and 



J3 = ^(12a+ + 56+) 

J3P3a = b+ (78) 
-'sPsfo = 0. 



When d+ = 0, 



When e+ = 0, 



Js = ^(12a+ + 56+ + 36e_ 

J3P3a = h+ + 12e+ 
J3P3h = 0. 



J3 = J(12a+ + 56+ + 144cZ+) 



(79) 



J3P3a = 6+ (80) 

J3p3h = 48(i+. 

If any one of a+, 6+, c+ or e+ is zero, we can determine a+, 6+, c+, d-^- 
and e+ from only the X3 distribution. However, if c?_|_ = 0, then we cannnot 
determine a+, 6+ and e+, independently. 

6.2 Results if either 7^ 0, 6+ 7^ 0, (i+ 7^ or e+ 7^ 
Similarly, When only a+ 7^ 0, 

J3 = 4a+ 

P3a = (81) 

P3b = 0. 



When only 6+ 7^ 0, 



P3a = I (82) 
P3b = 0. 
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When only a+, c+, (i+ / but a+ is negligible, 



J3 = 48d+ 

/03a = (83) 
P3b = 1- 

When only c+, e+ 7^ but 6+ is negligible, 

J3 = 12e+ 

JaPaa = 1 (84) 

-^3P36 = 0. 

In these case, psa and p^i, are independent from a+, 6+, d+ and e+. 

The results in subsection 14.21 and here is summed in the table below. 



Table 1: values of J' s and p's in some cases 





Ji 


Pla 


Plb 


h 


P3a 


P3fe 


(A) 


-ja+ 


3 
5 





4a+ 








(B) 


IK 


3 
14 


27 
28 




3 
5 





(C) 
















(D) 


24d+ 


1 





48(i+ 





1 


(E) 


18e+ 


1 

3 


1 


12e+ 


1 






(A) : all but gi, g2, gs, 94 are vanishing 

(B) : all but 55, g^ are vanishing 

(C) : all but eAji, eAi are vanishing 

(D) : all but Re[gseA*j^] + Re[g4eA*p^], cAr, eA^ are vanishing 

(E) : all but Relg^eA*^] + Re[gQeA*j^], eAji, eA^ are vanishing 



6.3 What happens ii t ^ jj^ dominates 

Similar to subsection 14. 3[ in this case, 



J3 = 12(4d+ + e+) 
Jspsa = 12e+ (85) 
JsPsb = 48(i+. 
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From Eqs. (f23ill and (f236]) in appendix El 



4 If J3P3b\ 

.2 , / 7. . N 2 

ia 



(86) 

+ - c+ c+ V 12 

So, we can determine a+ and 6+ lower limit, though we can't determine a+ 
and 5+, directly. 

We note here that since the relation Jspsf, = 48^+ is always true, a+ 
lower limit is always determined. 

6.4 Lower limit of c+ when r ^ is highly suppressed 

We give three types of c+ lower limits. Especially, these are very useful 
when r — > /i7 is highly suppressed compared with r ^ 3/i and c+ cannot 
determine directly. 

First, we give the relations, 



3J3 - 5 Jspsa - SJspsb = 12(a+ - 2e+) 
3J3 - 3J3/93a - SJspsb = 2(6a+ + b+) 
hPza = h+ + 12e+ 

from Eqs. (f69l) . 

Using Eqs. ([HTD and (pM]) and the fact that a+, h+> 0, 



(87) 



576( J3 - Js/Jsa - JsPsb) 6a+ + 6+ a 
Using Eqs. ([871) . 



192(3J3 - SJspsa - 2,Jzp3h) a+ - 2eH 



(89) 



Here, if c+ < 6+, then Jgpsa = 5+ + 12e+ = 0(6+) and 6+ > 0(1126+1). So, 
if (3J3 - 5J3P3a - 3J3P36)/12 > J3P3a, this means a+ - 2e+ > 0(6+) and 
then a+ — 2e+ ~ a+. So, using the relation ()234p . 

^-^''''^^ < c+. (90) 



192(3J3 - 5J3P3a - 3J3P3b) 
If 

3J3 - 5J3P3a - 3J3P3fe = 12(a+ - 2e+) < 0, (91) 



23 



6+ ^ c_|- and there is no special cancelation between a+ and 2e+, then 
a+ - 2e+ = C'(-2e+) and a+ < 0(2e+) < 6+/6. So, 

(3J3 - 5J3/53a - 3J3P35)2 _ (a+ - 2e+)2 _ ^ /4\ ^ ^^^^^^ ^g^^ 



864(J3 - J3/03a - -hp^b) 4(6a+ + 6+) 

We note here that 6a+ + 6+ is similar to hj^ + 12e+ in this situation. So, 
we have similar result if the denominator of left hand side in Eq. ([92]) , J3 — 
JzP^a - JsPsb = 2(6a+ + 6+)/3 is exchanged by IJ^p^a/S = 2(5+ + 12e+)/3. 

From ()77p and ()234p . if we suppose that = from some special models 
or other experimental results, 

1 iJsPsbf _ 4 < (93) 



242 J3 - Jspb 



Similarly, from (f76]) and ()236p , if we suppose that a+ = from some special 
models or other experimental results, 

1 (3J3-5W_4<,^. (94) 



864 J3 - Jspia b 

6.5 A bound for e+ 

6.5.1 Without using information c+ 

We obtain a bound for e+ using information from Gq{xi) measurement 
without using information on c+. 
Similarly to the subsection 14.61 



^P^- 12 ^+ .Q.N 



3 - 3/93a - 3/>3b 6a+ + 6+ 6a+ + 6+ 
and if e+ = 0, 

< ^ , < 1. (96) 

3 - 3p3a - 3p3b 

This becomes equation when 6+ = and = respectively. So, if 

P3a 



1 - /53a - P3b 



< 0, (97) 



then e+ < 0; if 



then e+ > 0; and if 



> I (98) 



1 - PSa - P3b 2 



< ^ < |, (99) 
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then, from Eqs. 



Furthermore, we can determine e+ allowed region as follows. Using the 

fact 

< „ \, < 1, (101) 

which become the equation when = and a+ = 0, respectively, 
12 2p3a ^ 12 



6a+ + 6+ 3 - 3p3a - 3p3b 6a+ + 6 
The solution of this inequality about e+ becomes 



3J3 + 5 J3P3a + SJsPsb ^ ^ ^ ■/3P3a (103) 



24 - ^ - 12 

This becomes equality, when 04. = and 64. = 0, respectively. 

6.5.2 Using information 

We obtain a bound for e+ using information from Go{xi) measurement with 
using information on c+. 

Similar to subsection 14.6.21 from Eq. (j69p . 



b+c+ + 12e+c+ = J3/03feC+. (104) 

This becomes 

J3P3b+c+ > el + I2e+c+ (105) 
since 6+c+ — > as explained in appendix [Dl Using c+ > 0, 

J3P3b > ^^ + gy_3g_ (106) 



C-L \ c 



The solution of this inequality about e+ becomes 



-6c+ - c+ y' ^ + 36 < e+ < -6c+ + c+ + 36 (107) 

This becomes equality when &+c+ — = 0. So we have another limit of e+. 
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7 Physics Implications from Gq{xi) and ^0(2:3) Dis- 
tributions 



From the analysis of both of Gq{xi) and Go (2^3) distributions, we can deter- 
mine not only a+ , 6+ , c+ , d+ and e+ but also in more detail in some suitable 
cases as explained in this section. 

Introducing new real parameters {^i, r2, rs, r4, rs, rg, r/j, r^} > and 
27r > {^1,^2; 03,04,05, (^6, 9r,Gl} > 0, the effective coupling constants can 
be explained as 



51 
93 
95 



rie 



105 



i0E 



92 
96 



r2e 



102 
104 



The 



Then, 



(r? + rl) 



1 

16 

2 I 2 



+ ri + ri 



-r^ri cos(6'3 - Ol) - r4r/jcos(6'4 - Oji) 
-rQVR cos(6'6 - Or) - r^VL cos{05 - Ol). 



7.1 What can we say if a+c+ — (i^ = 

From Eq. ()234p in appendix [Dl a+c+ — = only if 

Im[gzeAl] + /m[(74eA|j] = 
, 51 = 52 = 0. 



In that case, 



r^me ^ = r^TLe 

from the relation (jllOp . This means 



(108) 



(109) 



(110) 
(111) 
(112) 

(113) 



O3-Ol = O^-0r, Oi-0R±27r. 

Then, from the relation (|111|) . 

r^VL sin(6'3 - Ol) + t^tr sin(6'4 - 9r) 
= r^rL sin(6'3 - Ol) + r^rR sin(6'3 - Ol) 
= (rsrL + r^rR) sin(6'3 - Ol) = 0. 



(114) 



(115) 
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(116) 



This means 

r^ri + r4TR 
and/or 
03 - = 0, ivr. 

Here, r^^ri + r/^rpi ^ because of the conditions (|235p and (jll2|) . Then 
B^-eL = OA-OR = 0, zbvr and 

d+ = ^{r^rL + r^rR). 

The sign of right hand side is minus 63 — 9l - 
ez-9L = 9i-eR = ±7r. 
In these cases, 



(117) 
and plus if 



a+ = 
c+ = 
d+ = 
rsrR 



rl + rl 



(118) 



One of these conditions is dependent on others. For example, one condition 
c+ = r|j + r\ can be expressed using other conditions as 



4 

a+ 



(^r-i+r4rfl)2 
jt- 



+rl 



^R + ri 



(119) 



4 

since a+c+ — = 0. So, there are only three independent conditions in 
(jllSp . Using one of r^, r4, vr and r^, we can express others. For instance, if 
we know rR from other experiments or some special models, we can represent 
other coupling constants as 



rL 



r3 




(120) 



Here, we note that if {^3, r4, r/j, r^,} / 0, 

|2 



\d^\ 



193 1 



1941 



153 1 



\eAR\^ 

1541 



(121) 



c+ leAi] IcArI 

from the relation (jllUp . 

If rs = 0, then = from the relations (f235]l . (fTTOl) and (fTT2]) . Simi- 
larly, r/j = if r4 = 0; r3 = if r^, = 0; and = ii vr = 0. 
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7.2 What if it so happens that a+c+ — > 

When a+C-|_ — (i^ > 0, at least, one of the relations 



/m[53^2] + Im[g^A\] / 
51/0 
52/0 



(122) 



is satisfied. 

First, if q^Ar - g^Ai / 0, then 53 / and/or Ar / Al- This 
leads parity violation which is defined in Appendix[El Next, if I'm[g^A*^ + 
/m[54^|j] / 0, then the relative phases of one of 93^2 and/or gA-^*R have 
nonzero values. This leads CP violation. If gi / or 52 / 0, the scalar 
current and/or pseudo scalar current exists. 

7.3 What can we say if 6+c_|_ — = 

From Eq. (j236|) in appendix [Dj 6_|-C_|- — = only if 




(123) 
(124) 



In that case. 



r->rRe 



(125) 



from the relation ()123p . This means 




(126) 



Then, from the relation (I124|) 



r^ri sm{95 - Ol) + r^rR sin(6'6 - Or) 
= r^ri sin(6'5 - 9l) + rQVR sin(^5 - Ol 
= {r^fL + TQrR) sin(6'5 - Ol) = 0. 



(127) 



This means 



r^TL + rQTR = 
and/or 

05-01 = 0, ±7T. 



(128) 



28 



Here, r^r^ + r^rji / because of the conditions (I237p . Then 9^ — 0l 
06 — Or = 0, ibvr and 



e+ = TinrL + r^TR). 

The sign of right hand side is minus \i 9^ — 9^ - 
e5-9L = 9Q-9ii = ±7r. 
In these cases, 



(129) 

9r = and plus if 



rf + r^ 
rl + rl 



e+ = T{r5rL + rern) 



(130) 



One of these conditions is dependent on others. For example, one condition 
c+ = + r\ can be expressed using other conditions as 



(131) 



b7 = ^+ 



since 



0. So, there are only three independent conditions in 



(jl30p . Using one of rs, rg, tr and r^, we can express others. For instance, if 
we know r/j from other experiments or some special models, we can represent 
other coupling constants as 



re 



' R 




(132) 



Here, we note that if {^5, rg, r^, r^,} / 0, 



1551 



155 1 



be I 



(133) 



from the relation (|123p . 

If rs = 0, then r^, = from the relations p37p and (I123p . Similarly, 
= if re = 0; rs = if = 0; and rg = if r/j = 0. 
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7.4 What if it so happens that 6+c+ — > 

When — > 0, at least, one of the relations 

55^/? - ge^L + 
lm\g^Al\ + lm{g^A\\ / 

is satisfied. 

First, if 55^/2 - flfe^L / 0, then 55 / and/or / A^. This 
leads parity violation which is defined in Appendix[El Next, if lm\gr^A\\ + 
Im\g%A*^ 7^ 0, then the relative phases of one of 95^2 aiid/or g^A*^ have 
nonzero values. This leads CP violation. 



7.5 What if one of A^,, (73, g^g^, g^ — 

We reveal here that even if a+c+—ci^ > or 6+c+ — > 0, we can determine 
the lower limit of l^sl, {g^l, {g^l, {gel, {bArI and {eAil in the case Al = 0, 
Aji = 0, Al = 0, Ar = 0, gs or 35 = and g^ or g^ = 0, respectively. If we 
restrict that A^ = from other observation for example r decay or 

some specific models. 



(P 

^ = r|cos2(e4-0^) < |54|2, 

c+ 

e2 

-t =ri cos2 (06- ^i?) < |56|'. 
So, the lower limit of {g^l, {g^l is determined. Similarly, if Ar = 0, 

^ = ricos2(03-eL)<b3p, 
e2 

^ = ricos2(05-^L) < I55I'; 



if 53 = 0, 



rl 



a+ ^2 _^ 



' 4 ^ 16 

if 55 = 0, 

„2 



and if 54 = 0, 



^3 16 



(135) 



(136) 



-4cos2(04 - 0^) < r|cos2(04 - Or) < jeAi^^; (I37) 



^ = r|cos2(06-0^)<|e^fl|2; (138) 

0-1- 



'2 2 

^3 cos2(03 - 0i) < rl cos2(03 - 0i) < le^Lp; (139) 
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if 96 = 0, 



-t = rlcos\6,-eL)<\eALf. (140) 

0+ 



The lower limits of {g^l, {gel, \g3\, {g^l, {bArI and {eA^l are determined in 
each case. 

7.6 What can we say about parity and CP symmetries 

If parity or charge symmetry exists, 

ri = r2 
ra = r4 
r5 = re 

O3 = ±04 (141) 

6*5 = ±6*6 

rR = rL 
Or = ±9l. 

So, 

a+c+ - 4 = !M + 4riri sm^O^ - 6l) (142) 
from Eq. (f23il) . and 

6+C+ - = Arlrl sin^O^ - Ol) (143) 

from Eq. (1236]) . If a+c+ - d\ ^ then |5i| = jc/al ^ and/or CP is 
violated. If 6+c+ — / 0, then CP is violated. 
If both of parity and CP symmetries exist, 

ri = r2 

rs = r4 

''s = '•e 

03 = 04 = 0, TT (144) 
05 = 0g = O,7r 
J-R = rL 

eR = eL = ^,TT. 



So, 



a^c+-di = ^ (145) 



from Eq. (f23il) . and 

6+C+ - = (146) 
from Eq. (l236l) . When a+c+ - / 0, |c/i| = I52I / 0. 
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7.7 What happens if Br{T fi'-f) times the fine structure 
constant is much smaller than Sr(r —>■ 3//) 

In the case c+ <C others, we can't determine c+ directly. However, using the 
relations (psil) and (1236]) . 

(147) 

b+ 

So, we can determine c+ lower limit. 

This method is very useful to predict r — > /i7 branching ratio in the case 
that r — > 3// is detected while t ^ fi^ has not been detected, yet. 



8 Angular Distribution, Gf{xi) 

As we shall see, angular distribution of the particle a in r — > ^'T-+a+anything 
will give us information about a_, b-, c_, d-, e_, /+ and g+. 

The general formula Q contains all the information. However, it is 
too complex to analyze it here. So, we integrate about d^l to simplify the 
formula. Furthermore, we introduce three formulae, which are the general 
formulae integrated about three kinds of azimuthal angles of the momentum 
ka as explained in Figs. [7119] to simplify the formula, respectively. Thanks 
to this prescription, we obtain three simple formulae which are convenient 
to analyze here. 



frame 2 

z K 




Figure 7: The polar coordinate of frame 2. The North pole is z direction. 
Consider the double differential distribution, where d^ka is replaced by 
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dkad COS dzd(pz' 



da 



dxidx2di}rdip dyad(j)azd cos 9a 



Br{T — > iivv)Br(T v^- -\- a + anything) ° 
Go(xi,X2)G?(y„)(2 + -i) 



- G|(xi,X2)G§(y„){(2 - ^)kazPiz - P^kyPiy + (l + ^)kaxPix} 

(148) 

where ya = 2Ea/mr. X)i X2)Pi2, ^iGI{xi,X2)Pix, Y.iGl{xi,X2)Piy 

can be obtained by considering three single differential distributions as fol- 
lows. 

Integrate d(l)az, then 
da 



dxidx2dCl-,-d'(p dyadcosOa 



= 2Tr Br {t — > pLvv)Br{T — i^t + a + anything) 



Go(xi,X2)G?(y„)(2 + ^) 
-Y^Gl{xuX2)G'i{ya){2-^)kazPi 



(149) 



And using the equation 

da 



dxidx2dQ,rd'4} dya 



d (3y 

A'KBr{T — > iivv)Br(T + a + anything) " 



xGo(xi,X2)G?(ya)(2 + -i), 



then 



da 



1 



da 



dxidx2dQ,rdi^ dyad cos 9az 2 dxidx2d^rd'ip dya 



-2TrBr{T iii'i')Br{T Vt + o- + anything) " 

: GiiVa) (2 - ^)kaz Glixi,X2)Piz. 



(150) 



(151) 
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Similarly, Taking the coordinates as Figs. [8]and[9l we have the equations 
da 1 da 



dxidx2dQrdip dt/adcosOax 2 dxidx2dQrdip dya 
= —27rBr[T ^vv)Br{T^ ^ z^,- + a + anything) 



(152) 



and 



da 



1 



da 



dxidx2dftrdil^ duadcosOay 2 dxidx2dQ.T-dip dya 



2iTBr{T ^vv)Br{T ^ u^- + a + anything' 

X G'i{ya)P%ayY.Gl{xi,X2)Piy, 



(153) 



respectively. 

frame 2 



frame 2 





Figure 8: The polar coordinate of Figure 9: The polar coordinate of 
frame 2. The North pole is x direc- frame 2. The North pole is y direc- 
tion, tion. 

From the above three equations, we pull out the quantities, 

^PixGl{xi,X2) 
i 

^PiyGl{Xl,X2) 

i 

Y,nzGl{xi,X2). 

i 

First, using these quantities, we analyze Gf(xi,X2) and G2{xi,X2 
temporary. G|,^(xi,j;2) which is the ps component of ^jPjG|(xi, 



(154) 



j con- 
X2) is 
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defined as 

Gp^{xi,X2) = P3 • ^PiG-(xi,X2) = P3 • PlGi(xi,a;2) +P3 • P2G|(xi,X2) 
i 

= hx X] PixGl {xi , X2) + Pay ^ PiyGl {xi , X2) + ^ PizGl {xi , X2), 

i i i 

(155) 

where 



P3 • Pi = 1 - 2- 

X1X3 

P3 • P2 = 1 - 2 . 

aJ2a;3 

Integrating over X2 and defining the parameters 

6_ = bsP - bel' 

c_ = le^Rp - le^Lp 

o?_ = -(i?e[53eAl^]-i?e[54e^a) 

e_ = -{Re[gQeA*ji] - Re\g^eAl]), 

Gp^{xi,X2) becomes 



(156) 



(157) 



^13(^1)= / dX2Gl^{xi,X2) 



l—x\ 

1 



4(2x1 - l){a_(2xi - 1)(5 - Axi) + 12d_(2a;i - 1) + 3e_(3 - 2xi)} 



+ 8c. 



-{2(2x1-1) 



(xf - 13x1 + 13) 
(l-^i) 



- 3(2x? - 2x1 + 1) log — —] + 24 log [2(1 - xi 

Xl 

+ 6_|(8x? - 32x1 + 23)(1 - 2xi) - 24(1 - xi)^ log [2(1 - xi)] } 

(158) 

To determine the value of coefficient c_ , we define the function 

i^5(xi) = ^^^^G|3(xi). (159) 
The value of F5 where xi = 1 is 

F^{x^)l^^^ = c.. (160) 
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So we can determine the value of c_ . 

Then, to determine the values of a_, 6_, d- and e_, we define another 
function subtracting the term of coefficient c_ , 

Feixi) = 6{Gl^{xi) - (c_ term)} , (161) 

where 

(c_ term) = 8c-\ 2(2xi - 1)^^^-^ ^ - 

L {l-xi) 

- 3{2xl - 2X1 + 1) log [^^] + 24 log [2(1 - xi)] } . 

(162) 

We'll find that the 4 aspects of this function lead to determination of 
parameters a_, 6_, d- and e_. The value of Fq where xi = 1 is 

Fea = Feixi)\^^^^ = 4a_ + 5_ + 48d_ + 12e_. (163) 

The gradient of Fq where xi = 1 leads 

Feb = \4-F6{xi) = 36_ + 32d_. (164) 

6 dxi xi=i 

The gradient of Fq where xi = 1/2 leads 

Fee = \^Fq{xi) ,=-b- + 8e_. (165) 

6 axi xi=^ 

The integration value of Fq leads 



1 



Fed = 6 J dxiFsixi) = 6a_ - 6_ + 48d_ + 24e_. 



(166) 



Prom previous four equations, we can determine the parameters a_, 6_, d_ 
and e_: 



a_ 


= ^(-lOFea + 3F66 - 9Fqc + SFea), 


(167) 


fe- 


= ^(-6F6a + SFgfe - 3F6e + 4^6^), 


(168) 


ci- 


= ^(ISFga - 7^66 + QFqc - l2Fed), 


(169) 


e_ 


= 4(-6F6„ + SFeb - Fee + 4^6^). 
Id 


(170) 



Next, we analyze G^{xi,X2) using the relation 

GUxi,X2) = • VPiG?(xi,X2). (171) 



Pi X P2 
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Integrating it over X2, we define the quantity as 

4/+ {(2x1 - l)(2x? - 2x1 - 1) + 3(xi - l)xi log[^^] } 



1 — X\ 

Xl 



3(1 -xi^ 



Xl 

+ g+[ (2x1 - 1) (4xf - 10x1 + 7) + 6(xi - 1)2 log [^^] } 

(172) 



where 



/+ = -{lm[g^eAl] + Im[gieA*R]) 
g+ = -(/mffifge^Jj] + Imlg^eA}}) . 



Here, we define the function 



^r(..)^G^(xi)|^. 



(173) 



(174) 



Thus we determine /+ and 5+ from 



Fiixi) =-8/+ + 2(7+ 
2:1=1 

^^7(2:1) ^ =35+. 



(175) 



So we can determine the values of a_, 6_, c_, e_, /+ and g'+ 
separately. Adding this result to the result of the previous section i.e. 
a+, 6+, c+, (i+ and e+, we can determine the values of |(7ip/16 + l^ap, 
I52IVI6 + b4|', b5p, |56|', iiebs^I], Re[gAA\], Re[g^Al], 

Re[gQA\], Imlg^Al] + Im[giAW and Im[g^Al] + Im[gQA*^], separately. 
Moreover, in some suitable cases, we can determine \gi\ and \g^\ ( also \g2\ 
and \gi\), separately. To separate these are the main purpose of following 
subsection. 

9 Physics Implication from G^xa) Distribution of 

In sections [3] and El we studied that a+ , 6+ , c+ , and e+ are determined 
by the energy distribution. In section [HI we also studied that a_, 6_, c_, d_, 
e_, /+ and g+ are determined by the energy distribution and the angular 
distribution. 
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Here, we study some special cases using these parameters. Using Eq. 
(jlOSp . each parameter defined in (jl57p and (jl73p are represented as 



Qj— 




2 2 
'3 '4 




b- 


2 2 
= ^5 - ^6 






C- 


2 2 

= rR-rL 






d- 


= -r^VL cos(6'3 - 


- + r4r/jcos(6'4 


-0r) 


e_ 


= -rer/j cos(6'6 - 


- ^'fl) + ?'5rL cos(6'5 


-0l) 


/+ 


= -r^VL sm{03 - 


■ Ol) - r4r/jsin(6'4 ■ 


-0r) 


9+ 


= -r&VR sin(6'6 - 


- Or) - rsrL sm{05 ■ 





(176) 



9.1 A sufficient condition for existence of scalar and/or pseudo 
scalar current 

We make three types of relations which may reveal the existence of scalar 
and/or pseudo scalar current in this and following two subsections. 

First, in this subsection, we give a sufficient condition for existence of gi 
and/or (72. Using the fact that 



and the relation 

2 d\ + d'^ a+c+ - a_c_ 
^+ + 2 2 

= r|r| + r|r|j + 2r3rLr4rR sin{93 - Ol) sin(6'4 - Or) - ""^'^"^ ^ 

2 2 

= -Ti^^^ - tI^^4^ + 2r3rLr4rKsin(e3 - ^l) sin(04 - ^i?), 
16 2 lb 2 



(178) 



we can make a relation 
16 



T^2^2 _|_ ^2™2 / / ^2 / ^2 



- 16 " ^''^''^ '■3r-4 sin(03 - Ol) sin(04 " + V 1^ + V + 



a+c+ — a_c_ 1 yT"^ 2 w 2 2T f2 ^J- + 



-V(4-a2)(4-c2)-/:| 



(179) 



If 

"^"^ ; - y{a\ - al){c\ - cl) -fl- > 0, (180) 
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then 



|5i|'|e^L|' + |52|'|e^R|' >0. (181) 
So, in this case, we can get the result that gi and/or 52 7^ 0. 

9.2 What can we say about gi, g2 and/or CP Violation 

We now give relations which are convenience for determining the existence 
of scalar and/or pseudo scalar currents and/or CP violation of Lagrangian. 
One of them is 

If the left hand side becomes zero and a++a_ / 0, then sin(^3— ^2,) = gi = 
and (a+ + a-)/2 = If/sp. On the other hand, if the left hand side is larger 
than zero, at least, {g^l sin(^3 — ^i,) / or (/i 7^ 0. {g^l sin(^?3 — 7^ means 
CP violation, gi means that the scalar and/or pseudo scalar current 
exists. Similarly, we make the relation, 

- i^^)"^^ = -OR)+"f^>0. (183) 

If the left hand side becomes zero and a+— a_ 7^ 0, then sin(^4— ^r) = 52 = 
and (a+ — a_)/2 = |34p. On the other hand, if the left hand side is larger 
than zero, at least, \g4 \ sin(04 — 9r) 7^ or 52 7^ 0. I54I sin(6'4 — Or) 7^ 
means CP violation. 52 7^ means the that scalar and/or pseudo scalar 
current exists. 



9.3 Existence of gi and/or g2, or the Values of Im[g3eA}] and 
Im[g4eA*ji] 

Finally, we have two more relations. First one is 



a+ + a_ c+ — c_ /<^+ + d- „ /a+ — a_ c+ + c_ — d- „ 
- ( n ) + V Ti 7i ( ^ 



= rL\jrl sin2(^3 _ ^i) + + rj sm'^{04 - Or) + ^ 
> rsrLl sin(6'3 -9l)\+ t^vrI sin(04 - 9r)\ 



It becomes an equality when 



91= 92 = 

Im[g3eAl]Im[g4eA*j^] > 0. 



(184) 



(185) 
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In that case, from the relation ()184p and the conditions ()185p . Imlg^eA'^] 
and Imlg^eA"^] are determined as 



Im[g4eA*ji] 



' ^ : : (186) 



1/4 

If /+ = 0, it means Imlg^ieA^] = Im[gieA*j^ = since Im[g'ieA\]Im[gieA*^ > 
0. 

We derive another relation from (|178p 
2Im[g3eAL]Im[g4eAji] = + 



+ l{\g,WAL\' + \g2\'\eARf) 



16 



> i {2/2 - (a+c+ - a_c_) + (4 + )} . 



2 

When the equation of (I184p is not an equality and 2/^ — (a+c+ — a„c_) + 
(4 + c?^) > , then we can derive, using (|185p . that 51 / and/or g2 / 0. 
Furthermore, if 2/^ — (a+c+ — o_c_ 

) + (4 + '^-) > ' then /m[53eA2] / 
and /m[(74e74Jj] 7^ since I'm[g3eA*j^]Im[g4eA*j^] > 0. In that case, CP 
symmetry is violated and the sign of Im[g3eA*j^] and Im[g4^eA*j^] is the same 
as that of — /+. 

9.4 Imlg^eA}] and ImlgQcA}^] 

Here, we give the method to determine Im[g^eA^] and I'm[gQeA'^]. 
Substituting the relations 

Imlg^eAlf = Issl'le^Lp - Reig^eAl]^ 
Im[geeA*jif = IggPle^/jp - Re[geeA%]'^ 

to 

Imig^eAlf = {g+ + ImigQeA*^]^ 

= I'm[gQeA*j^f + 9^ + 2g+Im[gQeA*j^], 
the imaginary part of geeA*j^ is represented only by the observables as 
-ImjgeeA*^]'^ - gj + Im[g^eAlf 
25+ 

IffsPleAflp - Relg^eA^f + gl - bsPleA^P + Relg^eAlf 

25+ 

6+c_ - 6_c+ - 2e+e_ + 251^ 

(190) 



(188) 



(189) 



Im[gQeAW 
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So, if 5+ / 0, Im[gQeA*j^] can be determined independently. Similarly, 
Im[g5eA*j^] is represented as 

/m[55eA* ]2 + gl- Im[geeA*j^]'^ 



Imlg^eAl] 



2^+ 

Iffsl'leAj,!^ - Re[g5Alf+gl - \ge\^\eAR\^ + Re[g(ieA*j^]^ 

25+ 

6_c+ - 6+c_ + 2e+e_ + 25(i 



45+ 

if 5+ 7^0. 

Even if 5+ = 0, 

c-i- + c_ 6+ — 6_ / e+ + e ^ ^ 

2 2 V 2~ 

and 

c+ — c_ 6+ + 6_ / e+ — e_ 



(191) 



/mbee^a^ (192) 



= /mbse^I]^. (193) 
So we can determine the absolute values of Imlg^eA"^] and Im[gQeA'^]. 

9.5 What can we say if one of Ar, Al, g^, = 

If we restrict that Ar = from other experiments or some specific models, 
then 

/+ = -r3rL sin(6'3 - Ol) 

d+ = -rsVL cos(03 - Ol) (194) 



c+ = rl. 



So, 



"+ J+ 2 I |2 
'^+ 

a+ + a- _ dl + fl ^rj^lgil (195) 
2 c+ 16 16 

c+ = |e^Lp. 

In this case, I53I and |5i| are determined independently. Similarly, if ^4^ = 0, 

"+ ^ ■/+ _ |„ |2 

C+ 

a+-Q- _ 4 + /+ ^ (196) 
2 c+ 16 

c+ = {eAnl^; 
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if 54 = 0, 



(4 + /i)-^ = b3'^ 



2 



c+ — c_ 16 



(4 + /i)-^ = ^ (197) 



and if 53 = 0, 



2 ~ 16 ^ 



Q+ + a- _ |giP 
2 ~ 16 ■ 



9.6 CP violation 

The master formula ([9]) contains a part, 

(1 + COS^ T] - kazPsz - ff' SVC? 7] KyP^y 



1 /\ Sill 2?y /\ /\ 
+ (l + — ) sin^ ry ^a^-Psx {kaxPsz + KzPc 



(199) 



3xJ 



7' 

In this part, each term is proportional to 

kaziPl X pi)^, kayiPl X pi)j^, kaxiPl X pi)^, kax{Pl X Pi)^ + /Ca^(pi X pi)^ 

(200) 

Supposing CPT theorem, time reversal is equivalent to CP transforma- 
tion. In time reversal, momentum flips their directions. In fact, 

ka -K, Pi -Pl, P2 -P2 (201) 

In this transformation, (j200p flip their signs. This means ()199p part in ([9]) 
flips its sign. This part proportions to or /+. These coefficients are 
the imaginary parts of the interactions. So, if they do not vanish, then CP 
symmetry is violated. 

9.7 Ai,y4ij > others 

If aBr{T — > fij) is large compared to 4-fermi sector of Br{T — > 3/x), here 
a is the fine structure constant, it may be difficult to determine \gi\ to l^el 



42 



directly from r 3/x. It may be still possible to get a bound for these 
coupling constant from interference effects as follows. 

\93\^ =4 >rl coshes -Ol) 



94\'^ = ri>ricos\04-eR) 



|55|' = ri >ri cos2 (^5- 0i) 



|56p = ri > ri cos\ee - On) = (205) 
9.8 < others 

If aBr{T — jJL'y) is small compared to Br(r — 3/x), we can still get some 
lower bounds from the amplitudes of r — > /X7 decay as follows. 

\eAL 



\eAL 



\eAR 



10 Concluding Remarks 

Neutrino oscillation suggests that, flavor quantum number isn't conserved 
not only in quark sector but also in neutrino sector. KM ansatz implies that 
if neutrino sector violate the flavor quantum number, charged lepton sector 
also do it. In the Standard Model, this violation is too small to determine 
in any designed futur experiments. However, in some new models, it is 
suggested that this violation is going to be determine. If r 3/Lt event is 



{d+ + d.f 
2(c+ - c_) 

(d+ - d-f 
2(c+ +c_) 

(e+ - e.f 
2(c+ - c_) 



(202) 



(203) 



(204) 



.rlrlcos\e^-eL)_{d+ + d-f 



rR> 



rR> 



, 2 

TB + n 




2(a+ + a_) 


rlrl cos^(^5 — 


Ol) 


(e+ - e_)2 


rl 




2(6+ + 6_) 


rlr^ cos^(^4 — 


Or) 


(d+ - d_)2 


li , „2 
16 ^ M 




2(a+ — o_) 


r|r|cos2(06 - 


Or) 


(e+ + e_)2 



2(6+ - 6_) 



(206) 



(207) 



(208) 



(209) 
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detected, there are many models which are suitable to the first experimental 
result. However, at least, all models without one model is not true. So, if 
r — > 3/i event is detected, our analysis must be necessary to figure out if one 
model is allowed or forbidden. 

We assumed only Lorentz and gauge invariance of the Lagrangian and 
locality of the action. So, if we cannot fit the data to the differential cross 
section, it means violation of Lorentz or gauge invariance. 

From energy distributions (|5iP/16+|(73p) + (|g2p/16+|(74p) 
I^Rp + \Al?, Re[gAA*j^] + Re[g-iAl] and Re[gQA\] + Re[g^Al] can be de- 
termined. Using the angular distribution of decaying products of r+ and 
r", we can determine (biP/16 + l^sP), (b2p/16 + |54p), I55I, bel, \Ar\, 
\AlI Re[giA\l Re[g^All Re[g&A\l Re[g^Al], ImlgiA*^^] + Im[g^Al] and 
Im[g^A*^ + I'm[g^A\\ independently. We can determine the argument of 
relative phases, arg[(74AJj] and arg[g'3^2]) if -^"^[54^1?,] and /m[(73AJ^] have 
same sign, and if there are no scalar and pseudo scalar interaction. Simi- 
larly, we can also determine the argument of relative phases, arg[(5reAJj] and 
arg[5574^], if gj^ is nonzero. 

Even if r — > /x7 process is suppressed by a factor of 100 or more than 
r 3/1 process, we may still estimate the branching ratio of r ^ ;U7 from 
these interference: 



Br(r — > 3/i interference) 
Br(r ^7) 



g3 + g4 + fl5 + g6 
Ab^ + Al 



» 1, 



(210) 



where a is the fine structure constant. In concrete terms, for instance, even 
if |eAip is too small to determine directly, we may be still able to get a 



lower bound (d+ + (i_) /{2(a+ — a-)} as explained in subsection 19.8 
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A Fiertz Transformation 

In Eq. (I2|), the terms 

{fR(7af3fJ-L){P'R<7"^fJ'L), ^^^^^ 
{TLCral3l^R){P'L(^'^'^m), 
{fLCrafSfJ-R) (/ii?<7°^AiL) , 
{TR(7al3^J'L){P'L(^"'^ ^J-r), 

which are naively assumed particularly don't appear. The reasons are as 
follows. 

First, using the equations 

{'IplRa^'' ii2L ) {ipSRCrfMulpiL ) 

1 - - (212) 

= -6(V'lRV'4L)(V'3ilV'2L) + ■^{ll^lR'y^''lp4.L)i-il^3R<7fiu'ip2L) 
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and 

i'4'lRCr'"'ip4L)iip3RCr^j,u'4'2L) 



1,T „„. . . (213) 



where the tpiL = (1 - 75)V'j/2, ipiR = (1 + 75)V'j/2, 'ipiL = + 75)/2, 
V^i-R = V'i(l ~75)/2 and -i/^j where i = {1, 2, 3, 4} are the Dirac spinors. Then, 



i'4'lRCr'"''ll^2L)itp3RCrfj,u1p4.L) " {iplRCr'"' 'll^4L)ilp3RCrfiu1p2L) 
= -4:{'^lR'>p4L){i'3R'4'2L) + 4('0iijV'2L) (V'3rV'4l) 



(214) 



and also exchanging L and R, 



ii'lLCr^''i'2R)i'4'3LCr^,u1p4.R) - {tplLCr'"' ■4'4R){lp3LCr^,u1p2R) 
= -4:{'lpiLtpm){i'3L'>p2R) + 4:{i!iLi'iR)ilp3Ltp2R)- 



(215) 



Next, 

{i>lRa^"'lp2L){i>3LCrfMi'i^m) = ^ C'i(V'lijri'i/'4i?)(V^3LriV'2L) = ^216) 

i 

where Fj = {1,75,0-'^''} and Q are some coefficients. 
Finally, 

{i'lLtp2R){'>l'3Rtp4.L) = -'^{i'lLl^tp4R){i'3LliJ.'4'2R)- (217) 

Then, 



{lplLi'2R){tl'3R'^4.L) - (V'1L^4r) (V'3i?V'2L) 
= 2{i!iLl^'lp2R){lp3L'yfii>'lR) - 2{i!iLl^1piR){ljj3L'yfi'4'2R) 



(218) 



and also 

(.1plR'4>2L){i^3Li^4R) - (^li?V'4L)(V'3LV'2i?) ^^^g) 
= 2('0ii^7''V2L)('03iJ7M'^4L) - 2(V'l/?7^V'4L)('03iJ7MV'2L)- 

So it is proved that Eq. ([2]) is the general form of 4-Fermi type interactions. 

B Parts of Master Formula 

B.l production cross section of r pair with the polarizations 

Defining as polarization vector in frame 3 and 2, respectively, the 
differential cross section for the process e~^e~ — > 7* — > (s'^)t^ (s^) in the 
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center of mass frame, framel is [I2j: 



1 + cos^ ?7 + — ^) + (1 + cos^ r? 



sin^ r\ 



T 



2 - /? sin r/ s+s^^ 



7 



sin 2?^ 

7 



(220) 



where a ~ 1/137 is the fine structure constant, /5 = |p:^|/£', 7 = (1 — 
^2-j-i/2 _ j^Ijyi^^ E is the energy of or e~ in the initial state and |p!^| 
is the absolute value of momentum of r"*". As described in Fig. [3l 77 is 
the angle between the momenta of e"*" in the initial state and . = 
+Pg-)^) where p'^+ and p'^^ are the momenta of and e~ in the initial 
state, respectively, is the solid angle for the r"*" momentum. We note that 
the quantities, -E, |p!^|, p'^± and are defined in frame 1. 

B.2 differential branching ratio for r decay 

The differential Branching ratio for the process u^- + a + anything in 

the rest frame of r~ is jT3] ; 



dBr{T {s ) ^ Ur + a + anything) 

d^a 

Br{T^ Vt + a + anything) 



(221) 



where Gf{ya) and ^2(1/0) are the functions of ya defined in each a. These 
are written in the table 1 of Ref. [T3]. Here, 



Aa = y" dyaylGliya), 



(222) 



2Ea 



(223) 



Ea is the energy of a, = ka/|ka| and ka is the momentum of the particle 
a. 

We note here that physical vector quantities which we treat in this pro- 
cess are only s~ and k^. The only scalar made by these vector quantities is 
s~ • ka So, we can explain the differential branching ratio, Eq. (j22ip by only 
two terms which are proportional to G\{iia) and s~ • kaG2(ya)) respectively. 
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B.3 narrow width approximation 

The narrow width approximation is 



1 



vr 



|p2 _ (m-ir/2)|2 ~ mT^^^ 
B.4 the total Branching ratio 



m 



where — ^1. 
m 



(224) 



In Eq. ([TT]) . if = mx(4s) which is the Upsilon 45' mass, the differential 
cross section in the center of mass frame, frame 1 is 



da 



dxidx2 



, 647ra2 / 2ml 
Cfo{xi,X2)— I IH ^ — 



'1 



mi 



X Br{T ^uv)Br{T' 
Then, the total cross section is 



m 



Ami;. 

2 

T(4s) / V "^T{4s) 

z/'r + a + anything) . 



(225) 



1 i-(f )^ 



a 



da 



-dxidxs 



i-_fa 

^ 2 



647ra^ 



1 + 



2mi 



Ami 



dx\dx^ ^'^T(4s) y '^T(4s) 

X Br{T jjLVv)Br{T~ + a + anything) 

\a+ + \h+ + Ad+ + 2e+ + ^ (24 log [^] - is) c+ 



m: 



T(4s) 



(226) 



where 



13 



20. 



since 5 = 2m^/m^. Here, we derive the total cross section in center of mass 
frame of initial electron and positron since the total cross section is Lorentz 
invariant under the boost for the beam direction. So, this expression is able 
to apply to the B Factory experiment. 



m 



^G|,/(1927r'^), the differential branching 



C G's 

Using the width r(r ^ ^vu) 
ratio concerning the polarization of r is 

dBr(r+(s+) ^ ^i/X2^3) _ 3 ^ . ^ r.^ + s+ • P r^fr. ro^ 

dx^dX2dnrdi; - 2^^'^^ ^ f^Uu}[Go{x,,X2}+S P.G,(xi,X2. 

(227) 

Go{xi,X2) and Gf{xi,X2) are as follows. 
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C.l G'o(a;i,X2) 



Go(xi, X2) =4a+(2 - xi - X2){xi + X2 - 1) + &+((! - xi)xi + (1 - X2)x2) 

(1 - Xi)(xf + (1 - + (1 - X2){xl + (1 - X2)2) 



(1-Xi)(l-X2) 

+ 16d+(xi + X2 - 1) + 4e+(2 - xi - X2) 



C.2 Gf(xi,X2) 

G'i(xi, X2) = - 4a_xi(xi + X2 - 1) + - xi) 

xi (2xi(l - xi) + 2(xi + X2 - 1) - 1) 



(228) 



+ 4c_ 
-4d_ 
- 2e_ 



(1-Xi)(l-X2) 

xi(xi + X2 - l)(-xi + xi(x2 - 2) + 2) (229) 

(1-Xi)(l-X2) 
Xi((2-X2)(l-Xi)^+X2(l-X2)^) 
(1-Xi)(l-X2) 

G|(xi,X2) = G?(x2,xi) (230) 

2xi(xi - X2)X2 



G3(xi,X2) = -2/+(xi + X2 - 1) - g+{xi + X2 - 2) 



(X1-1)(X2-1) 

(231) 



D derivations of relations, a+c+ — (f^ and 6+c+ — 

Here, we introduce a convenient formulae a+c+ — and 6+c+ — to use 
in sections HI El El and O Using this formula, we can determine more about 
current structure. 
Using the relations 

«+c+ - \g3eAl + g4eA*jif 
= ^(r? + rl){rl + r|) + (r| + r|)(ri + r%) 

- i^lrl - rlr\ - 2r3rLr4rR cos^Os -61-0^ + Or) 
= rlr\ - Irsrjir^rL cos{93 + 9r-9a,- Ol) + rjrl + j^irf + rl){r\ + r|) 
= Vze^'^TRe^'^ - r4e^^Vie^^^|2 + ^(r? + vDirl + r|) 



\g3eAR - g^eAil^ + ^(rf + ri)(ri + r|) 



(232) 
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and 

\gzeAl + gieA\\^ -d\ 

= 4rl + r|r| + 2r3rLr4rR cos(^3 -OL-O4 + Or) 

- rfrf, cos2(6'3 - Ol) - Ir^rLrATR cos{03 - Ol) cosiOi - Or) - rjr^ cos2(6l4 - Or) 
= rlrl sin2(^3 - Bl) + ^r^rLVirR sin(03 - Ol) sin(04 - Or) + r|r| sin2(^4 - Or) 



(233) 



{rarL sin(6'3 - 0l) + rirR sm{0i - 6r)}^ 
[imig^eAl] + Im[gieA*ji])\ 



we give a useful relation 

a+c+ — 

= ^{rl + rl){rl + r|) + \ne''^rRe''^ - r^^'^L^'^f 

+ (rarL sin(03 - ^l) + r4rK sin(04 -Or)^ (234) 

= + |32p)(|eAL|2 + |eylij|2) + 1^36^^ - g^eAj:^ 

+ (/m[<73eA2] + Imyg^eA*^)" > 0. 

In this formula, there are three terms. Each of them has zero or positive 
value. First term (l^ip + |5'2p)(|e^Lp + |e^Rp)/16 has the information 
about scalar and pseudo scalar currents. Second term \g3eAR — g4eAL\'^ has 
information about parity symmetry. Third term (^Imlg^eA*]^] + Im[g4eA*j^])^ 
has information about CP symmetry. 

If a+ = 0, then ri = r2 = = = and also (i+ = 0. Similarly, if 
c+ = 0, then vr = vl = and also d+ = 0. In these cases, we can't use 
interference effect. So, we consider only the case 

when we use this formula. 
Similarly, 

b+c+ - el 

= Ir^e'^'TRe'^" - ree^^V^e^^^p + {r^rLsinie^ - ^l) + rer^ sin(^6 - ^i?)}' 
= \g5eAR - gQeAL\^ + [imlg^eAl] + Im[gQeA*R])^ > 0. 

(236) 

In this formula, there are two terms. Each of them has zero or positive 
value. First term \g^eAR — g^eAi^^ has information about parity symme- 
try. Second term [lm[g^eA*i\ + Imlg^eA*^'^ has information about CP 
symmetry. 
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If = 0, then ^5 = rg = and also e+ = 0. Similarly, if c+ = 0, then 
rji = tl = and also e+ = 0. In these cases, we can't use interference effect. 
So, we consider only the case 



(237) 



when we use this formula. 



E C, P, T and CP Transformation 

We define charge (C), parity (P), time reversal (T) and CP transformations. 
In Lagrangian, C transformation is 



95 ^ 
Ar ^ 



P transformation is 



91 ^ 
93 ^ 
95 ^ 
Ar ^ 
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94 

9l 
Al 
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94 
96 
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94 ^ 

9e ^ 
Al^ 



9i ^ 
93 ^ 
95- 
A*R- 



^ 9i 
-^95 
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94 

9l 
Al 
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(239) 



and T and CP transformations are 



91 
93 
95 
Ar 



9i 

* 

93 

* 

95 
A* 



92 
94 
96 

Al 



92 
94 
96 

A% 



(240) 



F Observables in e^e Center of Mass Frame 

The relation 

da f da 



daida2---dan J d/3i(i/32 • • • , 



1=1 



a,- — a,, 



(/3i,/32,--- ,[3m))df3idf}2-- 
(241) 



is useful for converting the variables from /?2, • ■ ■ ) /?m to ai, 02, • ' " ^otr, 
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To write it briefly, we define sr], crj, s(p, c(j), s9, c9, sip, cip, s^, c^, sx 
and cx as sinr/, cost/, sincj), cos(/>, sm6, cos 9, simp, cosip, sin^, cos^, sinx 
and cosx, respectively. Also we define the quantities p[,p'2,p'^, k'^ as the 
quantities pi,P2,P3, ka in frame 1. First, 
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Here, 6 and (j) which are defined in Fig. U] represent the direction of p3 in 
frame 3. r/ represents the direction of r"*" in frame 1. 
Next, 
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(243) 



where ?/' is the angle between P2-P3 plane and P3-2; plane as represented in 
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Fig. [5l Similarly, 
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Here, as Eq. (fTKHIl . 



cos ^ = 1 - 2 



1 - X2 



P3 • Pi 
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sin^ = \/l — cos^ ^ 



and 



1 — a^i 

cos X = 1 - 2 = P3 • P2 



X2X2, 



(246) 



smx 



C0S2 X, 



where sinx < since vr < x ^ 27r. 

Finally, we use the 3 ways to explain k'^ for the benefit of the simplic- 
ity of the phase space integral. The coordinates are as Figs. El [8] and [9l 
respectively. 



k„ 







CT] 





/ 1 




\ —sr] 

7 





( 





1 






—sr\ 



\ / 

sri 



cr? y 






7 

1 

1 

.-7/3 

7/5 \ / 



1 




sr/7 


cry7 / 



■7/5 \ 



7 

a I sin Oaz COS( 

;a| sin6la^ sine 

|ka| COS Oaz 



( \ 

|ka| sin cos <i)az 
|ka| sin Oaz sin (/)a^ 
I ka I COS ) 

\ 



(247) 





7 








-7/5 \ 






Ea 




\ 




— srj^f) 


CT] 





srj'y 






\ka\ COS 9 ax 














1 







|ka 


sin 0(12^ cos ( 


/'ax 




v 


—crj^P 


—srj 





c?77 / 


v 


\K 


sin 9 ax sin (; 








7 








-7/5 \ 


/ 








\ 






CT] 





sr]^ 




\K 


sin sin c; 


bay 












1 









ka cos 






v 




—srj 





cry7 / 


V 


\K 


sin 9 ay cos ( 


Pay 





(248) 



(249) 



53 



